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Abstract. We prove new velocity averaging results for second-order multidimensional equa- 
tions of the general form, C{'^x,v)f{x, v) = g{x, v) where C{\7x,v) := a{v) ■ \7 ^ — Vj • b('i;)V2;. 
These results quantify the Sobolev regularity of the averages, f{x, v)4>{v)dv, in terms of 
the non-degeneracy of the set {v : < 5} and the mere integrability of the data, 

(/, g) g {LP. J,, ^). Velocity averaging is then used to study the regularizing effect in quasi- 
linear second-order equations, C{VxiP)p — S{p) using their underlying kinetic formulations, 
xtv)xp = 9s- particular, we improve previous regularity statements for nonlinear con- 
servation laws, and we derive completely new regularity results for convection-diffusion and 
elliptic equations driven by degenerate, non-isotropic diffusion. 
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1. Introduction 

We study the regularity of solutions to multidimensional quasilinear scalar equations of the 
form 

j=0 j,k=l ■> 

where p : M is the unknown field and Aj, Bjk, S are given functions from M to M. 

This class of equations governs time- dependent solutions, p{t = xo,x), of nonlinear conser- 
vation laws where Aq{p) = p and Bj^ = 0, time-dependent solutions of degenerate diffusion 
and convection-diffusion equations where {-B^-fc} > 0, and spatial solutions, p{x), of degenerate 
elliptic equations where Aj = 0. The notion of solution should be interpreted here in an appro- 
priate weak sense, since we focus our attention on the degenerate diffusion case, which is too 
weak to enforce the smoothness required for a notion of a strong solution. Instead, a common 
feature of such problems is the (limited) regularity of their solutions, which is dictated by the 
nonlinearity of the governing equations. A prototype example is provided by discontinuous 
solutions of nonlinear conservation laws. In |LPT94a] . Lions, Perthame & Tadmor have shown 
that entropy solutions of such laws admit a regularizing effect of a fractional order, dictated 
by the order of non-degeneracy of the equations. In this paper we extend this result for the 
general class of second-order equations (jl.lj) . In particular, we improve the |LPl'94a] -regularitv 
statement for nonlinear conservation laws, and we derive completely new regularity results for 
convection-diffusion and elliptic equations driven by degenerate, non-isotropic diffusion. 

The derivation of these regularity results employs a kinetic formulation of (ll.lj) . To describe 
this formulation let us proceed formally, seeking an equation which governs the indicator func- 
tion, Xpix)iv) := sgn(f )(|/9| — \v\)+ associated with p, and which depends on an auxiliary velocity 
variable v e M, borrowing the terminology from the classical kinetic framework. To this end, 
we consider the distribution g = g{x,v), defined via its velocity derivative d^g using the formula 

(1.2) d,g{x,v) := (a(t;) ■ V, - Vj • b(y)V, + S{v)dJ)xpix){v), a,- := A'-, b.^ := B'^, > 0. 

Observe that the nonlinear quantities $(p) can be expressed as the f-moments of Xp? ^{p) = 
^'{v)xp{v)dv, $(0) = 0. Therefore, by velocity averaging of ()1.2|) we recover ()1.1|) . Moreover, 
for a proper notion of weak solution p, one augments ()1.1|) with additional conditions on the 
behavior of $(p) for a large enough family of entropies $'s. These additional entropy conditions 
imply that g is in fact a positive distribution, g = m E A4~^, measuring the entropy dissipation 
of the nonlinear equation. We arrive at the kinetic formulation of (jl.l|) 

(1.3) C{Vx,v)f{x,v) = d^m - S{v)dyXp(x)iv), f = Xp, m e 

where £ is identified with the linear symbol C{i^,v) := a{v) ■ i^ + (b(t>)^,^). We recall 
that p{x) itself can be recovered by velocity averaging of f{x,v) = xpi'^): via the identity 
p{x) = f f{x,v) dv. In Section 121 we discuss the regularity gained by such velocity averaging. 
There is a relatively short yet intense history of such regularity results, commonly known as 
'velocity averaging lemmata'. We mention the early works of [GLPS88], |DLM91j and their 
applications, in the context of nonlinear conservation laws, in |LPT94a^ lLPT94'Bl ILPS96| IJ P02j : 
a detailed list of references can be found in |Pe02j and is revisited in Section |21 below. Almost 
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all previous averaging results dealing with (jl.Hj) were restricted to first-order transport equa- 
tion, deg£(z^, •) = 1. In Section f2.2l we present an extension to general symbols £'s which 
satisfy the so called truncation property; luckily — as shown in Section EiH below, all £'s with 
deg>C(i^, ■) < 2 satisfy this property. If the problem is non- degenerate in the sense that* 

3«G(0,1), /3>0 s.t. snp\{v: \Cit^,v)\ < 6}\ < (^gY , 

then by velocity averaging of the kinetic solution / in p.3p . p{x) = f{x,v)dv has a ly*'^- 
regularity of order s < Pa/ {3a + 2) with an appropriate r = > 1; consult ()2.25|) below. In 
the particular case of non-degenerate homogeneous symbols of order k, 

3a G (0,1) s.t. sup|{t;: \C{i^,v)\ < 5} \ < 5" , 

l€l=i 

we may take [3 = k and velocity averaging implies a M/j^^-regularity exponent s < ka/{3a + 2). 
The main results are summarized in averaging lemmata 12.21 and 12.31 

In Section El we turn to the first application of these averaging results in the context of 
nonlinear conservation laws, pt + V- A{p) = 0, subject to L°°-initial data po- If the equation is 
non-degenerate of order a in the sense that 

\{v:\T + a{v)-^\<6}\<6'' and sup |a'(t;) ■ ^| < 5^"", V + 1^1^ = 1, 

{v. |T+a(D)-^|<<5} 

then for t > 0, the entropy solution p{t,x) gains Sobolev regularity p{t, ■) G W-^^^(M.'^) of order 
s < Sa = a/(2a + l). This improves the Sobolev-regularity exponent of order a/{a + 2) derived 
at [LPT94aj (while facing the same barrier of si = 1/3 discussed in |DLW05j ) . 

Section m is devoted to convection-diffusion equations. We begin, in Section I^TTl with second- 
order degenerate diffusion pt — '^d^^^^Bjk^p) = 0. The emphasis here is on non-isotropic 
diffusion, beyond the prototype case of the porous medium equation (which corresponds to 
the case when Bjk is a scalar multiple of the identity, Bjk = B6jk)- The regularizing effect is 
determined by the smallest non-zero eigenvalue A(f ) = A(b(t>)) ^0 of b(t>) := B'{v), so that 

(1.4) \{v:0<X{v)<6}\<6'' and sup |(b'(t;)^, 5^"", V |^| = 1. 

{v. A(t))<<5} 

Staring with initial conditions po G L°°, then the corresponding kinetic solution p(t > 0, •) gains 
Wj^'jl-regularity of order s < 2a/ {2a + 1). In Section we take into account the additional 
effect of nonlinear convection. The resulting convection-diffusion equations, coupling degenerate 
and possibly non-isotropic diffusion with non-convex convection governing capillarity effects, 
are found in a variety of applications. Consider the prototype one-dimensional case 

Pt + A{p), - B{pU = 0, A'{p) ~ B'{p) ~ IpI" . 

The regularizing effect is dictated by the strength of the degenerate diffusion vs. the convective 
degeneracy, li n < i we find PF^^'^^-regularity of order s < 2/{n + 2), which is the same Sobolev- 
regularity exponent we find with the 'purely diffusive' porous- medium equation, i.e., when 
^ = 0. On the other hand, if the diffusion is too weak so that n > 2i, we then conclude with a 
Sobolev-regularity exponent of order s < l/{i+2), which is dominated by the convective part of 
the equation. In Section we present similar results for multidimensional convection-diffusion 

*We use X ^ Y to denote the estimate X < CY where C is a constant which can depend on exponents such 
as a,l3,p and on symbols such as £, Aj^ Bjk but is independent of fields such as p, co-ordinates such as x,t,v, 
and scale parameters such as 5. We use X ^Y to denote the assertion that X <Y < X . 
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equations with increasing degree of degeneracy. In particular, consider the two-dimensional 
equation 

Pt + (9., + d,,)A{p) - - d.,,fB{p) = 0. 

If we set A = 0, the equation has a strong, rank-one parabolic degeneracy with no regularizing 
effect coming from its purely diffusion part, since (b(w)^,^) = 0, V^i — ^2 = 0, indicating 
the persistence of steady oscillations, pq{x + y); moreover, if i? = then the equation has no 
regularization coming from its purely convection part, since a(f ) • ^ = 0, V^i +^2 = indicates 
the persistence of steady oscillations pQ{x — y). Nevertheless, the combined convection- diffusion 
with A{p) ~ p^"*"^ and B{p) ~ |p|"p does have l^^'^-regularizing effect of order s < 6/ (2 + 2n — £) 
for n > 2i, consult coroUarv 14.51 below. 

Finally, in Section El we consider degenerate elliptic equations, —^^dl.^^Bjki^p) = S{p). As- 
suming that the non- degeneracy condition (jl.4p holds, then the kinetic formulation of bounded 
solutions for such equations, p G W^'^lD), have interior regularity of order s < min(a, 2a/{2a + 
1)). We conclude by noting that it is possible to adapt our arguments in more general setups, 
for example, when a suitable source term S{p) it added to the time dependent problem, when 
dealing with degenerate temporal finxes, dtAolp), > 0, or when lower-order convective terms, 
Va; ■ A{p), are added in the elliptic case. 



2. The averaging lemma 

We are concerned with the regularity of averages of solutions for differential equations of the 
form 

(2.1) £(V.,t;)/ = A^9f^?, A.:=(-A^)V2. 

Here, / = f{x, v) G Wl^^^{R^ x M^) and g = g{x, v) G L^g^(Rf x M^) are real-valued functions of 
the spatial variables x = {xi, . . . ,Xd) G M*^ and an additional parameter t> G M, called velocity 
by analogy with the kinetic framework, and £(V^,f) is a differential operator on of order 
< k, whose coefficients are smooth functions of v. 

The velocity averaging lemma asserts that if C{-,v) is nondegenerate in the sense that its 
null set is sufficiently small — to be made precise below, then the w-moments of f{x, ■), 

7(x) := [ f{x,v)4>{v)dv, 0gCo~, 

are smoother than the usual regularity associated with the data of /(x, ■) and g{x, ■). That 
is, by averaging over the so-called microscopic f-variable, there is a gain of regularity in the 
macroscopic x- variables. There is a relatively short yet intense history of such results, motivated 
by kinetic models such as Boltzmann, Vlasov, radiative transfer and similar equations where 
the f- moments of / represent macroscopic quantities of interest. We refer to the early works of 
Agoshkov Ag84| and Golse, Lions, Perthame and Sentis, pPSSSl IGLPS88] treating first-order 



transport operators with /, g integrability of order p = q > 1 . The work of DiPerna, Lions and 
Meyer, |DLM91j . provided the first treatment of the general case p ^ q, followed by Bezard, 
JBe94' , their optimality in |Li95j and an optimal Besov regularity result of DeVore and Petrova, 

iDypoQj- 

Extensions to more general streaming operators were treated by DiPerna and Lions in 
[DLSQa', IDLSQbj with applications to Boltzmann and Vlasov-Maxwell equations, and by Li- 
ons, Perthame and Tadmor in |LPT94a] with applications to nonlinear conservation laws 
and related parabolic equations. Gerard, |Ge9Qj together with Golse |iGG92j provided an 
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L^-treatment of general differential operators. A different line of extensions consists of ve- 
locity lemmata which take into account different orders of integrability in x and v, leading to 
sharper velocity regularization results sought in various application. We refer to the results 
in |DLM91j for / e Lp^(R^, LP^iR^)), g G Li^(R^, L'^^iRi)) and to Jabin and Vega |JVn8] for 
/ G W^^'P'{R^,LP^{Ri)),g G W^^^^'{R^,L'^^{R^)). Westdickenberg, |Wen2j . analyzed a general 
case of the form / G B^^ p^{Ri,U^{Ry)),g G B^^^^^iR^, U^{R^)). Jabin, Perthame and Vega, 
HEQS], \TVQ^ used a mixed integrability of / G LP(Ri,W^^^P{R^)), g G Li(Ri,W ^^'''(R^)) to 
improve the regularizing results for nonlinear conservation laws [ LPT94aj . |LPT94b] . Ginzburg- 
Landau, and other nonlinear models. Golse and Saint-Raymond, |GSR02j have shown that a 
minimal requirement of equi-integrability, say / G L^(M^, L*(M^)), G L^{Ri x R^) measured 
in Orlicz space L* with super-linear is sufficient for relative compactness of the averages /, 
which otherwise might fail for mere L^-integrabihty, |(;LPS88j . 

The derivation of velocity averaging in the above works was accomplished by various methods. 
The main approach, which we use below, is based on decomposition in Fourier space, carefully 
tracking f{C,,v) in the "elliptic" region where {v : C{iC,,v) 7^ 0}, and the complement region 
which is made sufficiently small by a non-degeneracy assumption. Other approaches include the 
use microlocal defect measures and H-measures |Ge90t ITaQOj . wavelet decomposition [ DVPOOj 
and "real-space methods" — in time |Va991 IBD99j . and in space-time using Radon transform 
[TMnllW^ . X-transform |Wn8irTVn^ and duality-based dispersion estimates |C;SR()2j . Al- 
most all of these results are devoted to the phenomena of velocity averaging in the context of 
transport equations, k = 1, 

Our study of velocity averaging applies to a large class of £'s, satisfying the so-called trun- 
cation property: in Section 12.41 we show that all £'s of order k < 2 satisfy this truncation 
property. In particular, we improve the regularity statement for first-order velocity averaging 
and extend the various velocity averaging results of the works above from first-order transport 
to general second-order transport-diffusion and elliptic equations. The results are summarized 
in the averaging lemmata 12.11 and 12.21 for homogeneous symbols and in averaging lemma 12.31 
for general, truncation-property-satisfying £'s. Our derivation is carried out in Fourier space 
using Littlewood-Paley decompositions of / G W[^^{Ri x R^) and g G L'^od^t ^ ^v)- To avoid 
an overload of indices, we leave for future work possible extensions for more general data with 
mixed (x, v)-integrability of / and g. 

2.1. The truncation property. We now come to a fundamental definition. 

Definition 2.1. Let m(^) be a complex-valued Fourier multiplier. We say that m has the 
truncation property if, for any locally-supported bump function on C and any 1 < p < 00, 
the multiplier with symbol i(j{m{^)/6) is an L^-multiplier uniformly in 5 > 0, that is, its LP- 
multiplier norm depends solely on the support and size of ip (for some large i which may 
depend on m) but otherwise is independent of 6. 

In Section 12.41 we will describe some examples of multipliers with this truncation property. 
Equipped with this notion of a truncation property, we turn to discuss the L^-size of parame- 
terized multipliers. Let 
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where f{^,v) = T^f{i,v) : 




/(x, v) dx 
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denote the operator associated with the complex-valued multiplier m{-,v), truncated at level 
S. Our derivation of the averaging lemmata below is based on the following straightforward 
estimate for such parameterized multipliers. 

Lemma 2.1 (Basic estimate). Let J be a finite interval, / C and assume m{^, v) satisfies 
the truncation property uniformly in v E I. Let 1 < p < 2. Let Ai^ denote the velocity- 
averaged Fourier multiplier 

M^f{x) := j^M^f{x,v) dv = j^J^~^ip{m{^,v)/5)J^^f{x,v) dv. 
For each ,^ e M'' and 5 > let Qm{C,;S) C / be the velocity set 

6):=[vel : ^^^^ e supp 

Then we have the multiplier estimate 

1 

(2.2) \\M^^f{x)\\LP(^Rd) < sup |fim(^; 5)]^ ■ ||/||LP(RdxR„) 

Proof. For p = 2, the claim (j2.2j) follows from Plancherel's theorem and Cauchy-Schwarz in- 
equality, while for p close to 1+ the claim follows (ignoring the bounded \ flm{^] 5)|-factor) from 
the assumption that m(^,f) satisfies the truncation property uniformly in v (and in fact, the 
end point p = I, with the usual Ti^-replacement of L^, can be treated by a refined argument 
along the lines of |DLM91j ). The general case of 1 < p < 2 follows by interpolation. □ 

Remark 2.1. Clearly, if m(^, ■) and c(^) are L^-multipliers then so is their product, and in 
particular, if m has the truncation property then ()2.2|) applies for ilj{m{^,v) /6)c{^). 

2.2. Averaging lemma for homogeneous symbols. We will present several versions of the 
averaging lemma. The later versions will supercede the former, but for pedagogical reasons 
we will start with the simpler case of homogeneous symbols. In this case it is convenient to 
use polar co-ordinates ^ = |^|^', where ^' G S'^~^ is defined for all non-zero frequencies ^ by 

We begin with 

Averaging Lemma 2.1. Let I < q < 2 and let g G Lj^^(Rf. x M.^) if q > 1 or let g be a locally 
bounded measure, g G M{Ri x M^) if q = I. Let r],N >0 and let f G W^^^{Ri x M„), a > 
0, 1 < p < 2 solves the equation 

(2.3) C{V,,v)f{x,v) = kld^g{x,v) in V'{Wi x M,). 

Here, C{Wx-,-) is a differential operator with sufficiently smooth coefficients, C{-,v) G 

and let C{i^,v) be the corresponding symbol. We assume that C{i^,v) is homogeneous in ^ of 

order k, k > a + rj, that the modified symbol C{i^\v) obeys the truncation property uniformly 

in V, and that it is nondegenerate in the sense that there exists an a, < a < (A^ + l)q' , such 

that 

(2.4) sup \nc{^-S)\<S'', nc{^;S):={veI : \m,v)\<s}. 
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Then, there exist 6 = 6a & (0, 1) and Sa = s{6a) > cr such that for all bump functions 
(f) G C^{I), the averages f{x) := / f{x,v)(f){v)dv belong to the Sobolev space Wf'^^iM.^) for 
all s G (cr, Sa) and the following estimate holds 

(2.5) ||7lkfj;(R^) ^ lklll^^^(RdxR,) ■ ll/ll^{;P(KdxK„)' seia,Sa), s« := (1 - + ^«(/c - r/). 
Here 9 = 6a{p, q, N) and r are given by 

^ ' a{l/p' -l/q') + N +1 r p q 

Remark 2.2. It would be more natural to assume that the symbol C{i^,v) itself, rather than 
the modified symbol C{i^',v), obeyed the truncation property, as it is typically easier to verify 
the truncation property for the unmodified symbol. Indeed when we turn to more advanced 
versions of the averaging lemma (which rely on Littlewood-Paley theory, and which do not 
assume homogeneity) we will work with the truncation property for the unmodified symbol. 
However we choose to work here with the modified symbol as it simplifies the argument slightly. 

Proof. We start with a smooth partition of unity, 1 = z) such that ipo is a bump 

function supported inside the disc 1^1 < 2 and the other tpj^s are bump functions supported on 
the annulus 1/2 < \z\ < 2 (we note in passing that the other ipj^s can be taken to be equal, so 
the index j merely serves to signal their 'action' on the shells, 2^^^ < \z\ < 2-'+^). We set 

/,(x, v) := V, {^^^§f^) m v), J = 0, 1, 2 ... , 

recalling that ^' := ^/|^|, and we consider the corresponding decomposition f = fo + J2j>i fj- 
We distinguish between two pieces, / = f^^^+f^^^ where f^^^ := /o and f^^^ := J2j>i fj- Observe 

that the w-support of /(°) is restricted to the degenerate set Qc{^;6) whereas /(^) = X]j>i /i 
offers a decomposition of the non-degenerate complement, ^^£(^; S). The free parameter 6 is to 
be chosen later. 

We start by noting that f^'^^ = fo is associated with the multiplier ipo(C{i^' , v) / S) . Since 
C{i^',v) satisfies the truncation property, we can use Lemma ITT] and non-degeneracy assump- 
tion ()2.4|) . to obtain 

(2.7) ||7W||^.,.(^,) < sup nci^;6) 



ICI=i 



We turn to the other averages, /j, j > 1 which make /"^^^ Since £(i^, ■) is homogeneous of 
order A;, equation ()2.3|) states that 

and thus, 

^l-'^h = ^^K' jj^i^^r-y^aii, v)<P{v)dv, J = 1, 2, . . . , 
where ipj{z) := il)j{z)/ z is a bump function much like ipj is. Integration by parts then yields 

(2.8) Ar'':^= m.v)mdv+ 

+ lower order or similar terms. 
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We can safely neglect the lower or similar order terms which involve (powers of) the bounded 
multipliers dlC{i^',v),i < N and we focus on the leading term in ()2.8|) . associated with the 
multipliers 



(2.9) ^<">(£(|^)£^"(,r,„). 

By Hormander-Mikhlin or Marcinkeiwicz multiplier theorems, {iC,', v) are bounded multipli- 
ers and hence ()2.8p are upper-bounded by 

1 

{2^6)^ 

~{N) 

Here Aij = Ai~(N) are the Fourier multipliers with symbol ipj {C{i^' ,v)/2^). 

We now fix g > 1. By our assumption, C{iC,', v) satisfies the truncation property, and Lemma 
12.11 implies (^d) < (2-5'+^5)"/'? \\g\\L<i (^d). Adding together all the f/s, we find that 

/(^) = Ej>i fj satisfies 

(2.10) WAt^miLl^im) ^ E 7^7^(2^'^'^)^ II^II^L«xM.) ^ SV-^''^'^\\9\\Ll^m.u.y 
Thus, if we fix t > and choose 5 to equilibrate the bounds in (j2.7j) and (j2.1Up . 

then this tells us that 



i nf 

/(O) +/(!)=/ 



and the desired PF^^'J-bound follows for s < (1 — 6)a + 6{k — rj) with 9 given in ()2.(ip . The 
remaining case of g = 1 can be converted into the previous situation using Sobolev embedding. 
In this case, g being a measure, it belongs to W~''''^'' for all (e, such that 

d + 2 ^ d + 2 

geW <e<l<q,< — — , 

q'^ d + l 

and hence ()2.5|) applies for s < {l — 6)a + 6{k — rj — e) and 6 = 9a{p, qe, N); we then let e approach 
0+ so that q^ approaches arbitrarily close to 1+ to recover ()2.5|) with Sa and 6a{p, 1, N). □ 

Remark 2.3. As an example consider a (possibly, pseudo-) differential operator C^Vx, ■) of order 
k and let f{x,v) G Wj^^^ such that C(Vx,v)f G Assume that Cy{i^,v) 7^ so that 

the nondegeneracy condition ()2.4|) holds with a = 1. Application of the averaging lemma ITT] 
with p = g = 2,A^ = and rj = k — a — 1 then yields the gain of half a derivative, f{x) G W[^^ 
with s < (j/2 + (cT +l)/2 = cT-l-1/2, in agreement with [GG92t Theorem 2.1]. The main aspect 
here is going beyond the L^-framework, while allowing for general and possibly different orders 
of integrability, {f,g) G {^^^,11). 



Remark 2.4. The limiting case of the interpolation estimate ()2.5|1 . 6 = I, s = Sa, corresponds 
to Besov regularity / G B^f^iM'^). This regularity can be worked out using a more precise 
bookkeeping of the Littlewood-Paley blocks. For the transport case, = 1, it was carried 
out first in |DLM9H Theorem 3], improved in |Be94j and a final refinement with a secondary 
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index t = p can be found in IDVPOOj . This limiting case is encountered in the particular 
situation when k = t], so that the interval {a, Sa) 'survives' at ^ = 1. Here, one cannot 
expect for regularizing effect, but there is a persistence of relative compactness of the mapping 

To gain a better insight into the last averaging lemma, we focus our attention on the case 
where / is a ly^^f-solution of 

(2.11) C{\/,,v)f{x,v) = d,g{x,v), g E LL(Mf x M,),/ e W.^^^i^'' 



. X ^v), 



corresponding to the special case r] = and = 1 in ()2.3p . This case will suffice to cover all 
the single- valued applications we have in mind for the discussion in Sections 01 Eland El without 
the burden of carrying out an excessive amount of indices. The averaging lemma ()2.1|) implies 
that f{x) has Sobolev regularity of order s < {1 — 9)(T + 9k, 

(2.12) Il7lk-{R^) ^ lklllf„^(R^xR„) ■ WfWluRiv o = e^:= ^^Y/^7ZT7^^T2- 

The last regularity statement can be improved. To this end, we revisit the dyadic multipliers 

^ -(AT) 

in ()2.9|) . ipj [C{i^' ,v)/2^^C^ {i^' ,v). The key observation is that C^iii'^v) acts only on the 
subset of f 's — those which belong to f G ^ci.^'] Linking the size of C^ii^' ,v) to that 

of C{iC,',v), we arrive at the following improved averaging regularity lemma. 

Averaging Lemma 2.2. Let f E W^jof (I^f xl^tj), o' > 0,1 < p < 2, solves the equation 

(2.13) C{y:„v)f[x,v)=d^g[x,v), ^ ^ j A<(M^ x M^), q = l. 

Let v) be the corresponding symbol. We assume that C{i^, v) is homogeneous in iC, of order 
k,k > a, that the modified symbol C{iC,',v) satisfies the truncation property uniformly in v, and 
that it is nondegenerate in the sense that there exists an a, < a < q' , such that Ii2.4\ ) holds. 
Moreover, assume that 

(2.14) 3/i G [0,1] s.i. sup sup \C^{i^,v)\ < S'', nci^;5) := {v E I : \C{i^,v)\ < 6} . 

Then, there exist 9 = 9a E (0, 1) and r given by 

(215) 9-= 1.= !^ + ^ 

^ ■ ^ ■ a(l/p'-l/g') + 2-/i' r ■ p q 

such that for all bump functions cj) E C^{I), the averages f{x) := J f{x,v)(f){v)dv belong to the 
Sobolev space W['^^{R^) for all s E (a, Sa), Sq, = (1 — 9a)cr + 9ak and holds. 

For the proof we revisit ()2.8|) with 77 = 0, = 1, 



Its bound in (j2.1Up can now be improved by the extra factor of (2^~^^S)^ which follows from 
yielding 

and we conclude by arguing along the lines of the averaging lemma I7!T1 □ 
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Remark 2.5. In the generic case of a symbol C{iC,, v) which is analytic in uniformly in v, the 
'degeneracy of order a' in ()2.4|) implies that sup u)| : |^| = 1, \C{i^,v)\ < 6} ^ 6^~°'. 

Thus, fl2.14|) holds with /i = 1 — a and in this case, the velocity averaging holds with Sobolev- 
regularity exponent 



f{^) e W^f^iK), s < := (1 - e^)a + e^k, 9^ := 



a{l/p' + l/q) + l' 



Remark 2.6. The statement of the averaging lemma IT^ can be extended for / G Lp(R'^ x M„) 
and g G ^^^^(R^ x M^) in the full range of 1 < p < cxd and 1 < g < oo. The dual claim to (|2.2|) . 
based on (L^, BMO) interpolation for 2 < p < cxo reads 



(2.16) \\M^f{x)\\Lp(^^d) < sup IVlmi^; 6)\P ■ ||/||LP(RdxR„), 2 < p < oo. 

This yields the same iy*''^-regularity as before 

(2-17) WJWw^YA^^) ~ ll^llLf^^(M^xIR„) ■ ll/llw|{f(RdxK,)' * ^ ~ ^°)^ + ^"^) 

with 9 = 9a (p, g) given by 

a/p 11-6'^ _ , n - f i\ 

9:=—— — , -:=^^ + =, p := max(p,p), g := max(g,g). 

q;(1/P — 1/g) + 2 — /i r p g 

2.3. An averaging lemma for general symbols. We now turn our attention to averages 
involving general, not necessarily homogeneous symbols; as such the polar co-ordinate repre- 
sentation ^ = |^|(^' is no longer useful and will be discarded. We focus on equations of the form 
£(Va;, f )/(a;, f ) = d^gix^v)^ corresponding to {ri,N) = (0, 1) in (|2.3|) . 

Averaging Lemma 2.3. Let f G VF^^f (Mf xM^), o">0,l<p<2, solves the equation 

(2.18) C{V:,,v)f{x,v) =d^g{x,v), ^ ^ | A^(Rf x M^), g = 1. 

Let C{i^, ■) be the corresponding symbol of degree < k with sufficiently smooth v-dependent 
coefficients and assume it obeys the truncation property. Denote 

u;c{J;S):= sup \nc{^;6)\, Qc{^;6) := {v. \C{z^,v)\<6} 

and suppose the following non- degeneracy condition holds 

(2.19) 3a,/5>0 s.t. ujciJ; 6) < (^-j^Y , V5 > 0, J > 1. 
Moreover, assume that 

(2.20) 3X>0 and [0,1] s.t. sup sup |/:^(z^, t;)| < J^V. 

Then, for all bump functions G C^{L), the average f{x) := J f{x,v)(j){v)dv belongs to the 
Sobolev space VF;^'J(Mf ) for s G (a, Sa,i3) and the following estimate holds 



(2.21) 



f{x, v)(f>{v) dv 



w,. 



||/(x,tO||vy-.p(MdxM.) + l|£/(a;,^^)||LL(RdxR„). 
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Here, Sa,i3 '■= (1 — Oa)<^ + 0al3{2 — fx — X) where 6 = 6a and r are given by 

a/p' 11-99 
2.22 0:=—-, '-f- , - := + -, < ^ < 1. 

^ ' a{l/p' -1/q') + 2- fi' r p q 

Remark 2.7. How does the last averaging lemma compare with the previous ones? We note 
that since deg£^(z^, ■) < fc, then the additional assumption ()2.20p always holds with /i = and 
/5A = k. Hence, the averaging lemma with just the non-degeneracy condition ()2.19p yields 
the regularity f{x) G Wif^{W!j) of (the reduced) order s G (a, s^,/?), where 

Oil v' 

(2.23) Sa,p := (1 - 9a)(y + 9^(2^ - k), 9^ := . , ^f, 

a{l/p' - l/q'j + 2 

Now, if in particular, C{i^, ■) is homogeneous of degree k, then 

-) r = ^ 

this shows that if the non-degeneracy condition (j2.4p of averaging lemma ITT) holds. uJc{J] ^) ~ 
((^/J'^')", then it implies (|2.19|) with f3 = k, and we recover the homogeneous averaging lemma 
I2.1l with [rj, N) = (0, 1), namely, the averages / gain regularity of order s < {l—9)a+9{2j3—k) = 
(1 — 9)a + 9k. The only difference is that now the truncation property is assumed on the 
unmodified symbol C{i^,v) rather than the modified one C{i^',v). 

As for the averaging lemma 12. 2t we first note that in the generic case of an homogeneous 
£(-,f), the additional assumption ()2.20p holds with \ = a and /i = 1 — a, 

(2.24) sup sup \C^{i^,v)\ < J^'^d'^-''. 

|?|~J {vel: \C{i^,v)\<S} 

Indeed, all the homogeneous examples discussed in SSections 13 141 and El below, employ the 
averaging lemma 12.21 with these parameters which yield VT^'^'-regularity of order s < (1 — 
9a)a + 9aP, 



(2.25) f{x)eW[f^{Rt), s<{l-9a)a + 9af3, = ^i/^, + i/^) + ^ 

In the particular case of C being homogeneous of order k then j3 = k and we recover the 
averaging lemma l272l (except that the truncation hypothesis is now assumed on the unmodified 
symbol). 

Proof. We begin by noting that we can safely replace the non- degeneracy condition ()2.19|) with 
a slightly weaker one, namely 

(2.26) 3a,/? > 0, and £ > s.t. uJc{J]5) < [-j^^ , V(5 > 0, J > 1, 

and still retain the same gained of regularity of order s < Sa^p- This can be achieved by 
replacing the values of a in ()2.19j) by a/{l + e) and then absorbing e into a slightly smaller 
order of regularity, (1 —e)sa,i3 dictated by 9a/{i+e)- The extra e-power of 5 will be needed below 
to insure simple summability, which probably could be eliminated by a more refined argument 
involving Besov spaces, along the lines of |DLM91j . 
Next, we break up / into Littlewood-Paley pieces, 

/ = /o + E 

dyadic J's^Jl 
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SO that the spatial Fourier transform of fj{x,v), is supported for frequencies |,^| ~ J, 

and /o has support in |^| < 1. Since /o is a smooth average of / at unit scales, the contribution 
of /o is easily seen to be acceptable. By giving up an e in the index Sa,f3 we may thus reduce 
fj2.21|) to a single value of J. It thus suffices to show that 



J'WfjWLl^iR") ~ ^'^ll/j||Lf„^{RrfxR„) + Wghl^iRixR,), fj ■= / fj{x,v)(f){v)dv, 

J V 

for each J > 1. 

Fix J. Because of the local nature of Littlewood-Paley projections when J > 1 we may 
replace the localized L'^ norms with global norms. Actually we may replace U' by weak U' 
since we may pay another e in the index Sa,i3 to improve this. By the duality of weak and 
U''^ it thus suffices to show that 

(2.27) \{Tj,Xe)\^J-'\E\'/^' 

for all sets C M'^ of finite measure, where we have normalized ||/j||LP(iRdxM„) < and 

We now decompose the action in v-space of each of the Littlewood-Paley pieces (rather than a 
decomposition / itself used in lemma 

fj{x,v)= 2^ )fA^M. ^\ )-=^x'l^\ — — 

dyadic &'s1^J^ 

Here, '^/'(-z) is a bump function on C supported on the region \z\ ~ 1. It suffices to estimate 



(2.28) ( / ^(^^^)iA-Mmdv.XE) 



with a summable decay as the dyadic 5^0, so that ()2.27|) holds. 

By our assumption, C{iC,,v) satisfies the truncation property uniformly in v, hence by Lemma 
12. II we see that for all 1 < p < 2, 

2.29) / ^(^i^l^]/,(a;,^)0(^;)rft; < o^^^; 5) II xR„) • 



5 

From ()2.29|) and Holder we may thus estimate ()2.28p by 



(2.30) |( y ^(^^^^^)/.(-,t;)0(t;)rft;,Xi.)| ^ J-ujc{J;6y/^'\E\'/^' . 
On the other hand, thanks to equation (j2.18p we can write 

where ip{z) := iIj{z)/z and the gj's are the corresponding Littlewood-Paley dyadic pieces of g. 
We thus have 

(2.31) I ^{^^^^)fj{x,v)<P{v)dv=^- I ^(^^^)|^^,(^,^)0(^)rf^- 

We now integrate by parts to move the d/dy derivative somewhere else. We will assume that 
the derivative hits il){C{W xiv) / 5), as the case when the derivative hits the bump function 0(f) 
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is much better. We are thus led to estimate 

(2.32) i|(y" ^^,(^i^^)£,(V.,t;)^?,(^,^)0(^)rf^,Xi.) 

Since gj is locahzed to frequencies ~ J, then by ()2.20|) . the multipher C^{i^,v) acts hke a 
constant of order 0{J^^6'^). Also, ipz is a bump function much like ip. Thus we may modify 
()2.29|) — with p replaced by q (assuming that g > 1 and using the modified argument for the 
case g = 1 as before), ip replaced by ipz, and fj replaced by Cy(Vx,v)gj, to estimate ()2.HH) by 

(2.33) |( I i^{^^^^)fA;v)<^iv)dv,XE)\ ^ 5-^''^^ J^'^J; 5)'^'' \E\'/'' ■ 
Interpolating this bound with ()2.30|) . we may bound ()2.28|) by 

The parameterization in ()2.22|1 dictates a/r' = 6(2 — /i). Finally, we put the extra e power into 
use: by fl2.26|) . coci-J] SY^^' ~ (5"'^'^^ J"^)^*"^ and hence the last quantity is bounded by 

^e(2-M)ej-s„,,|^|l/r'^ = (1 - 9)a - 9^(2 - /i - A). 

Summing in 6 and using the hypothesis that s < Sa^p we obtain ()2.27p . □ 

2.4. Velocity averaging for first and second order symbols. To apply the velocity aver- 
aging|231"we need to find out which multipliers m(^) have the truncation property. Fortunately, 
there are a large classes of such multipliers. 

First of all, it is clear that the multipliers m(^) = ^ ■ ei and m(^) = have the truncation 
property, as in these cases the Fourier multipliers are just convolutions with finite measures. 
Now, observe that if m(^) has the truncation property, then so does m(L(^)) for any invertible 
linear transformation L on M'^, with a bound which is uniform in L. This is because the L'^ 
multiplier class is invariant under linear transformations. 

Because of this, we see that the multipliers 

mi(0 = a{v) ■ 

and 

have the truncation property uniformly in f , where a(f ) are arbitrary real coefficients, and b(v) 
is an arbitrary elliptic bilinear form with real coefficients. 

From the Hormander-Mikhlin or Marcinkeiwicz multiplier theorems and the linear transfor- 
mation argument one can also show that mi(^') has the truncation property uniformly in v. 
These arguments go back to the discussion of |DLM91j . The situation with m2(^') is less clear, 
but fortunately we will not need to verify that these second order modified multipliers obey 
the truncation property since our averaging lemmata also work with a truncation property 
hypothesis on the unmodified multiplier. 

Now we observe that if mi(^), m2{C,) are real multipliers with the truncation property, then 
the complex multiplier mi{^) +im2{^) also has the truncation property. The basic observation 
is that one can use Fourier series to write any symbol of the form 

'mi(0 +im2(0' 
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as 

where ipj{x) := e^'^^^^ip and ip is some bump function which equals 1 on the one-dimensional 
projections of the support of ip. Since the norm of tpj grows polynomially in j and ipi^j, k) 
decays rapidly in j, k (if ip is sufficiently smooth), we are done since the product of two 
multipliers is still an U' multiplier. 

3. Nonlinear hyperbolic conservation laws 

Having developed our averaging lemmata, we now present some applications to nonlinear 
PDE. We begin with the study (real- valued) solutions p(t, x) = p{t, xi, . . . , Xd) € L°°((0, oo) x 
K^) of multidimensional scalar conservation laws 



d 




in P'((0,oo) X Rf). 



We abbreviate ()3.1|) as + ■ A(p) =0 where A is the vector of C^'^-spatial fluxes, A : = 
{Ai,A2,...,Ad). 

Let X7(f ) denote the velocity indicator function 

( 1 ifO<t;<7 
X^{v) = I if 7 < < 

I otherwise. 

We say that p{t, x) is a kinetic solution of the conservation law (j3.ip if the corresponding 
distribution function, Xp{t,x){'v), satisfies the transport equation 

(3.2) dtXp{t,x)i^) + a(t;) ■ V.Xp(m-)(^) = dMt,^,^) m ^^'((0, oo) x x R,), 

for some nonnegative measure, m{t,x,v) G 7W + ((0, oo) x R^ x R^,). Here, a(t>) is the vector 
of transport velocities, a(t>) := (ai(f ),..., arf(f)) where aj(-) := A^(-),j = 1,2, ...d. The 
regularizing effect associated with the proper notion of nonlinearity of the conservation law p.ip 
was explored in |LPT94a] through the averaging properties of an underlying kinetic formulation. 
For completeness, we include here a brief description which will serve our discussion on nonlinear 
parabolic and elliptic equations in the next sections and we refer to |LPT94a] for a complete 
discussion. 

The starting point are the entropy inequalities associated with (|3.1|) . 

dMp{t,x)) + V,-A\p{t,x))<0 in P'((0,oo)xRf). 

Here, rj is an arbitrary entropy function (a convex function from R to R) and A'' := {A^, . . . , A^) 
is the corresponding vector of entropy fluxes, A^{p) := ri'{s)A'j{s)ds, j = 1,2, . . . d. A 
function p G L°° is an entropy solution if it satisfles the entropy inequalities for all pairs 
{rj, A"^) induced by convex entropies r]. Entropy solutions are precisely those solutions which are 
realizable as vanishing viscosity limit solutions and are uniquely determined by their L°° fl L^- 
initial data, po{x), prescribed at t = 0, e.g., |Lax73j . A decisive role is played by the one- 
parameter family of Kruzkov entropy pairs, {rj{p] v), A^{p; v)), parameterized by f G R, 

V{p; ^) := \p-v\, A]{p; v) := sgn(p - v){Aj{p) - Aj{v)). 
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Kruzkov entropy pairs lead to a complete L^-theory of existence, uniqueness and stability of 
first-order quasilinear conservation laws, |Kr70j . We turn to the kinetic formulation. We define 
the distribution, m(t,x,v) = mp(^t,x){v) by the formula 

(3.3) Mt. ^. := - k"^"' - + V, ■ (^"'^^ - 

The entropy inequalities tell us that the distribution m = nip is in fact a nonnegative measure, 
m(t,x,v) G 7Vl+((0,oo) X X R^). Next, we differentiate ()3.3|) with respect to v: a straight- 
forward computation yields that Xp{t,x){v) satisfies the kinetic transport equation ()3.2|) . This 
reveals the interplay between Kruzkov entropy inequalities and the underlying kinetic formu- 
lation; for nonlinear conservation laws, kinetic solutions coincide with the entropy solutions, 
jPT911 ILPT94aj . Observe that by velocity averaging we recover the macroscopic quantities 
associated with the entropy solution p, 

Xpiv)(l)iv)dv = <l>(p) 

where := J^^q (j){s)ds is the primitive of 0. In particular one can recover p itself by setting 
4>{v) = 1[_M,M](^^), M = \\p\\l°°- 

We now use the averaging lemma ()2.2j) to study the regularity of p. To this end we first 
extend (j3.2j) over the full x M^-space, using a C^(0, oo)-cut-off function, = 1 for t > e, so 
that f{t,x,v) := Xp{t,x){v)^{t) and g{t,x,v) := mp{t, x,v)'tp{t) + Xp{t,x){v)dt^{t) satisfy 

dtf{t,x,v) + Si{v)-VJXt,x,v) = d,g{t,x,v) in P'(Mt x x R,), (7 G ^(R^ x x M,). 

Set / := [inf po, sup po] and assume that the first-order symbol is non-degenerate, ()2.4j) . ()2.14|1 . 
namely, 

(3.4) 3aG(0,l)s.t. sup |na(C; 5)1^ 5", n^{^;6) := {v e I : \t + ai{v) ■ ^\ < d], 

r2 + 15|2 = l I. J 

and 

(3.5) 3p G [0, 1] s.t. sup sup \a'{v) ■ ^\ < 5^. 

151=1 !^a(€;<5) 

We apply the averaging lemma 17!^ for first-order symbols, k = 1, with g = l,p = 2 and 
a = 0, to find that J{x) and hence p{x) belong to VT/^c (^t ^ 

Pit, ^) e WZ{{e^ 00) X ), . < := " r := " + ^ ~ 

a + 4: — zp a + Z — p 

At this stage, we invoke the monotonicity property of entropy solutions, which implies that 
for s < 1, ||p(t, ■) II w^'Vr'*') nonincreasing, and we deduce that p(t,-) G W^'^{R'^) for t > e. 
We conclude that the entropy solution operator associated with the nonlinear conservation law 
(EH), (El), Po(-) ^ p{t, ■), has a regularizing effect, mapping L°°(Mf ) into iy';](Rf), 



Vt > e > : p{t, ■) G Wi'^^iX), s < Si, Si := 9. 



a + A-2p 

Next, we use the bootstrap argument of |LPT94al §3] to deduce an improved regularizing 
effect. The W';'^ (Rt X R^)-regularity of pit, x)'4){t) implies that fit, x, v) = Xp{t,x) {'v)'^{t) belongs 
to L^{W'''^{Rt X moreover, since d^Xpiv) is a bounded measure, / G L^iW^'^iRt x 

Mf),M^) nLi(Rt X M^, Vr^'i(R^)) and hence 
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Interpolation with the obvious L°°-bound of / then yields that / G Wi^^^(Rt ) for all 

s < Si/2. Therefore, the averaging lemma EiH applies to / = Xp{t,x){v)ip(t) with q = l,p = 2 and 
a = Si/2, implying that p{t, •) has improved VF^^'J-regularity of order s < S2 = {l — 6a)si/2 + 6a- 
Reiterating this argument yields the fixed point t Sqo = 2^^/(1 + ^a) and we conclude with 
a regularizing effect 



Vt > e > : poEL^n L\Ri) ^ pit, ■) G W^f^iRt), s < 



a 



a + 2 — p 

As indicated earlier in remark IT7| in the generic case, p = 1 — a, 
(3.6) sup sup \ai{v) ■ ^\ < 5^"°, 

151=1 {t,g/;|r+a(i.)-5|<5} 

which yields 1^'^' ^-regularizing effect of order s < a /(2a + 1). This improves the previous 
regularity result |LPT94a( Theorem 4] of order s < a/{a + 2), corresponding to p = 0. We 
can extend the last statement for general L^^^ initial data. Recall that the entropy solution 
operator associated with p.lj) is L^-contractive. We now invoke a general nonlinear inter- 
polation argument of J.-L. Lions, e.g., |TaOOt Interpolation Theory, Lecture 8]; namely, if a 
possibly nonlinear T is Lipschitz on X with a Lipschitz constant Lx and maps boundedly 
Yi I— > Y2 with a bound By, then one verifies that the corresponding i^-functionals satisfy 
K{Tx,t; X,Y2) < LxK^XjiBy/ Lx] X,Yi), and hence T maps [X, Yij^^g > [X,Y2]0^q. Conse- 
quently, the entropy solution operator maps [L^, L^]e,q ^ [L^, W^ioc]e,q, < 6* < 1 < g, and we 
conclude 

Corollary 3.1. Consider the nonlinear conservation law subject to U' fl L^-initial data, 

p{0,x) = po{x). Assume the non- degeneracy condition of order a, JJ.^I ), i3. 6|) holds over 
arbitrary finite intervals I. Then p{t,x) gains a regularity of order s/p', 



Vt > e > : poeL^n L^R^) ^ p{t, ■) G lyri(Mf), s < 



{2a + l)p'' 



The study of regularizing effects in one- and two-dimensional nonlinear conservation laws has 
been studied by a variety of different approaches; an incomplete list of references includes 

We close this section with three examples. Let i > 1 and consider the one-dimensional 
conservation law 

(3.7) |''^^'^) + ^{7Tl'''^'^^''^^}=°' /^o^C-^'^]- 

It satisfies the non-degeneracy condition ()3.4|) with a = hence p{t,-)\t>e £ ^lol with 
s < a /{2a -|- 1) = !/(£ -|- 2). It is well-known, however, that the entropy solution operator of 
the inviscid Burgers' equation corresponding to £ = 1, maps L°° BV , |U163j . This shows 
that the regularizing effect of order a /{2a + 1) stated in corollarv 13.11 is not sharp (although 
the averaging argument is! consult [DLWOSj following jJF02j). Accordingly, it was conjectured 
in |LPT94aj that (j3.4|) yields a regularizing effect of order a. 

Next, let m > 1 and consider the two-dimensional conservation law 
d , d ( I ,. .^ d 



dt 



VELOCITY AVERAGING AND REGULARITY OF QUASILINEAR PDES 17 

If £ 7^ m then ()3.4j) is satisfied with a = min{j, ^} and we conclude p{t > e, ■) G ^^^^(L^) with 
s < min{^, ■:;^^}- If i = m, however, then there is no regularizing effect since t'+v^^[+v"^^'2 = 
for r = 0,^1 + ^2 = 0; indeed, po{x — y) are steady solutions which allow oscillations to persist 
along X — y = const. Other cases can be worked out based on their polynomial degeneracy; for 
example, 

+ — sin(p(t,x)) + — |-p3(t,a;)| =0, po E [-M,M], 
has a non-degeneracy of order a = 1/4, yielding VT^^'J-regularity of order s < 1/6. 

4. NONLINEAR DEGENERATE PARABOLIC EQUATIONS 

We are concerned with second-order, possibly degenerate parabolic equations in conservative 
form 



(4.1) _p(t,x) + A,(p(t,x)) - d^^Mt.^)) = 0, in P'((0,oo) x M^). 

j=l ^ j,k=l ^ ^ 

We abbreviate, pf +Vx- A(p)+trace(V2:(8>V2: B(p)) = where B is the matrix B := 

By degenerate parabolicity we mean that the matrix B'(-) is non-negative, (B'(-),^, ^) > 0, V.^ G 

W^. Our starting point are the entropy inequalities associated with ()4.1|) . such that for a// convex 

(4.2) 9tr7(p(t,a;)) + V,-A^(p(t,x))-trace(v,®V, B''(p(t,x))) <0 in D'((0, oo) x ). 

Here, A'' is the same vector of hyperbolic entropy fluxes we had before. A'' = (A^, . . . , A^) and 
B'' is the matrix of parabolic entropy fluxes, B'' := (-Bjfc)j_fc=i5 -^Jfc(p) •= f v'{^)Bjki^)d'S. We 
turn to the kinetic formulation. Utilizing the Kruzkov entropies, ri{p;v) := \p — v\, we deflne 
the distribution, m(t,x,v) = mp{t^x)iv), 

■ Ti{p-v)-r^{Q-v) fA^{p- v) - A^(0; v] 

Ot r y X 



m(t,x,v) := — 

(4.3) + trace (^V, ® ' ^ 

The entropy inequalities tell us that m{t,x,v) E 7Vl+((0, oo) X X M^) and differentiation 
with respect to v yields the kinetic formulation. 



(4.4) dtXp{t,x){v) + a{v) ■ VxXpit,x){v) - • h{v)VxXpit,x){v) = d^m{t,x,v), 

for some nonnegative m E Ai^ which measures entropy+dissipation production. Here, a 
is the same vector of velocities we had before, a = A', and b is the non-negative diffusion 
matrix, b := B' > 0. The representation ri{p) — ri{0) = J T]'{s)xp{s)ds shows that the kinetic 
formulation ()4.4|) is in fact the equivalent dual statement of the entropy inequalities ()4.2|) . But 
neither of these statements settles the question of uniqueness, except for certain special cases, 
such as the isotropic diffusion, Bjk{p) = B{p)Sjk > 0, e.g., |(ya99j . or special cases with mild 
singularities, e.g., a porous- media type one-point degeneracy, [DiB93l ria97j . The extension 
of Kruzkov theory to the present context of general parabolic equations with possibly non- 
isotropic diffusion was completed only recently in |CP03j . after the pioneering work ^VH69j . 
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Observe that the entropy production measure, m, consists of contributions from the hyperbohc 
entropy dissipation and the parabohc dissipation of the equation, m = mA + fTi-B- The solutions 
sought by Chen and Perthame in |CP03j . p G require that their corresponding distribution 
function Xp satisfies ()4.4|) with a restricted form of parabohc defect measure me: the restriction 
imposed on me reflects a certain renormalization property of the mixed derivatives of p (or more 
precisely, the primitive of ^yh{p)). Accordingly, we can refer to these Chen- Perthame solutions 
as renormalized solutions with a kinetic formulation (j4.4|) . These renormalized kinetic solutions 
admit an equivalent interpretation as entropy solutions, |CP03j and as dissipative solutions, 
|PS05j . A general L^-theory of existence, uniqueness and stability can be found in |(yK04j . For 
a recent overview with a more complete list of references on such convection-diffusion equations 
in divergence form we refer to |Ch03j . The regularizing effect of such equations, however, is less 
understood. In |LPT94a( §5] we used the kinetic formulation ()4.4|) to prove that the solution 
operator, po ^ p(t, ") is relatively compact under a generic non- degeneracy condition 

sup 0)1 = 0, nc{^;0):={v : r + a{v) ■ ^ = 0, {h{v)^,O = 0}. 

ICI=i 

A general compactness result in this direction can be can be found |Ge90j . We turn to quantify 
the regularizing effect associated with such kinetic solutions. We emphasize that our regularity 
results are based on the 'generic' kinetic formulation ()4.4|) . but otherwise, they are independent 
of the additional information on the renormalized Chen-Perthame solutions encoded in their 
entropy production measure m. The extra restrictions of the latter will likely to yield even 
better regularity results than those stated below. We divide our discussion into two stages, 
in order to highlight different aspects of degenerate diffusion, in Section 14.11 and the coupling 
with nonlinear convection, in Section 

4.1. Non-isotropic degenerate diffusion. We consider the parabolic equation 



(4.5) -p{t, x) - J2 0^,^^'^^^^' "^^^ = °' 

j,k=l 

Here, we ignore the hyperbolic part and focus on the effect of non-isotropic diffusion. The 
corresponding kinetic formulation (j4.4j) extended to the full x x reads 

dJit, X, t;) - Vj • hiv)VJ{t, X, v) = d,g{t, x,v), f := Xp^(t), 9 G ^+(Ri x x M,). 

Set / := [inf po, sup po]- The corresponding symbol is £(r, i^, v) = it + (b(t>)^, and it suffices 
to make the non-degeneracy assumption ()2.4|) . on the second-order homogeneous part of the 
symbol C{Q,iS,,v) = (b(t>)^,^). We make 



(4.6) 3a e (0,1) s.t. sup |r]b(e;5)U ny,{^;S) := {v e I : < {h{v)^,0 < 5}. 

ISI=i 

and 

(4.7) 3pG[0,l] s.t. sup sup \{h'{v)^,0 \ ^ 5^" 

151=1 



We apply the averaging result IT^ with q = l,p = 2,a = and A; = 2, to find that f{t, x) and 

C(RtxMf), 



hence p{t,x) belong to Wi^dRt x ^"^^ 
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We follow the hyperbolic arguments. The kinetic solution operator associated with ()4.H1 is 
L-'^-contractive, hence ■)IIw'''Vr'*) nonincreasing and we conclude that Vt > e, pit, ■) has 

loc ^ X ' 

W;^'^ -regularity of order s < si, si := 26^ = 2a/ {a + 4 — 2p). We then bootstrap. Since 
Xpiv)'ip{t) G W[^^{Rt X X Ry) for all s < si/2, we can apply the averaging lemma 12.11 
with a = si/2 leading to VF^^'J-regularity of order S2 := (1 — 6a)si/2 + 26*0 with fixed point 
Sk T =4^^/(1 + 0,), 

(4.8) Vt>e>0: p, e n 1^^^) ^ p{t, ■) eW^^'M), s< . 

a + 2 — /i 

We distinguish between two different types of degenerate parabolicity, summarized in the 
following two corollaries. 

Corollary 4.1. [Degenerate parabolicity I. The case of a full rank] Consider the degenerate 
parabolic equation subject to L°° fl -initial data, p(0,x) = po- Let Xiiy) > \2{y) > 

...\d{v) > be the eigenvalues of h{v) and assume that Xd{v) ^0 over I = [inf po, sup po] . 
Then, {|b(t>)^,^) > Ad(f)|^P and p{t,x) has a regularizing of order s < 2a/{a + 2 — p)p' , i.e., 
holds with a and p dictated by Xd{v) = Xd(h{v)), 

|nA(^)U5° and sup sup | (b'(t;)^, ^ U (5^, nx{6) := {v e I : 0<Xd{v)<6}. 

\^\=iven^{5) 

Corollary 14.11 applies to the special case of isotropic diffusion, 

(4.9) ^p(t, x) - AB{p{t, x)) = 0, B'{v) > 0, 

subject to L°° n L^-initial data, p(0, ■) = po. If &(■) := B'{v) is degenerate of order a in the 
sense that, 

|fife(5)|<r, and sup |6'(t;)| < 5^"", Qb{6) := {v e I : < b{v) < 6} , 

then corollary mU implies Vt > e : p(t, ■) G W'^^, s < 2a/ {2a + 1). For Lp n L^-data po, the 
corresponding solution p{t, ■) gains M/^^'^^-regularity of order s < 2a/ (2a -|- l)p' and we conjec- 
ture, in analogy with the hyperbolic case, that the non-degeneracy (j4.6|) yields an improved 
regularizing effect of order 2a/p' . Existence, uniqueness and regularizing effects of the isotropic 
equation (j4.9p were studied earlier in fBCTD^ .BCSla^ .BC81b] . The prototype is provided by 
the porous media equation, 

(4.10) _p(t,a;)-A{-^|p"(t,x)|p(t,x)} =0, p(0, x) = po(x) > 0, po G 

The velocity averaging yields W^^'^-regularity of order 2/(n -|- 2) and, as in the hyperbolic case, 
it does not recover the optimal Holder continuity in this case, e.g. |DiB93| ITa96j . In fact, 
the kinetic arguments do not yield continuity. Instead, our main contribution here is to the 
non-isotropic case where we conjecture the same gain of regularity driven by Xd{h{v)), as the 
isotropic regularity driven by b{v). 

We continue with the more subtle case where b(-) does not have a full-rank, so that 

3i, l<i<d : Xi{v) >... Xe{v) > 0, Xi+,{v) = ... = Xd{v) = 

Despite this stronger degeneracy, there is still some regularity that can be 'saved'. To demon- 
strate our point, we consider the 2D case. 
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Corollary 4.2. [Degenerate parabolicity II. The case of a partial rank] We consider the 2D 
degenerate equation 



(4.11) 



d_ 
dt 



p{t,x) 



02 02 q2 

—Bu{pit,x)) + ^^-^B^Mt.x)) + —BMt.x)) 



0. 



dxi dxidx2 

subject to L°° nL^-initial data, p{0,x) = Pq. Assume strong degeneracy, bf2{v) = 4:bii{v)b22{v) , 

so t/iat A2(b(i;)) = 0, Vf G / = [inf po, supo] . In this case, {h{v)^,^) = (^Vbii{v)^i + \/b22iv)^2^ 

and p{t, ■) admits a W[^l -regularity of of order s < 2a/{a + 2 — p) which is dictated by the non- 
degeneracy. 



(4.12) 

and 
(4.13) 



sup 

l?l=i 



sup sup 

1^1=1 venucs) 



b'22iv) 

\/hf2{v 



-6 



We distinguish between two extreme scenarios. 

(i) If |6ii(f)| ^ 1^22(^)1 Vf G /, then the regularizing effect ()4.8|1 holds with (a,/i) dictated 
by b22{v), 

1^622(5)! < r and sup 1622(^^)1-^", Vti,.^M ■= [v ^ I ■ ^ < h22{v) < 5] . 



(ii) If 611 (v) = b22{v) Vt; G /, then there is no regularizing effect since the symbol y\/bii{v)^i + 

\/522('y)^2 ) vanishes for all ^1 ± ^2 = (so that ()4.12j) is fulfilled with a = 0). Indeed, the 
equation 1)4.111) . with bii{v) = 622 ("w) =: B'{y), takes the form 



and we observe that Pq{x =f?/) are steady solutions which allow for oscillations to persist along 
X ^ y = const. 



d' 



4.2. Convection-diffusion equations. We begin with the one-dimensional case 
(4.14) 



^p(t,x) + ^A(p(t,x)) - ^5(p(t,x)) = 0. 



We consider the prototype example of high-order Burgers' type nonlinearity, a{v) := v^,i > 1 
combined with porous medium diffusion b{v) = |u|",n > 1. The corresponding symbol is given 
by C{{t, i^),v) = iT+v^i^ + We study the regularity of this convection-diffusion equation 

using the averaging lemma 12.31 which employs the size of the set 

VLc{J]5) := J\T + v^i\ + .P\v['i'^ r^ + e^ = 1, J> 1, (5 < 1. 

Comparing diffusion vs. nonlinear convection effects, we can distinguish here between three 
different cases. Clearly, nc{J]S) C n^, := {v : 1^1" < 6/J^}, hence uJciJ;S) < (5/j2)i/" and 
()2.19p holds with at, = 1/n and /J^ = 2. We shall use this bound whenever n < i, which is the 
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case dominated by the parabolic part of ()4.14|) . Indeed, in this case we have (5/ jy^^ ^ (5/ J^)^/" 
which in turn yields 

sup \C,{{t,i^),v) \ < sup (J|i;|^"^ + J>r-^) < 
venc{J;5) veUb 

This shows that ()2.20p holds with {Xb,fib) = (ofo, 1 — and velocity averaging implies Wi^^- 
regularizing effect with Sobolev exponent s < (3hab / {"io^b + 2), 

2 

S < Sn = • 

2n + 3 

We also have, VLc{J] 5) C fi^ := {v : \v\^ < S/J}, so that uJc{J] 5) < {5/jy/\ i.e., fICTIl holds 
with aa = 1/i and f3a = ^- We shall use this bound whenever n > 2i, which is the case driven 
by the hyperbolic part ()4.14|) . In this case, ^ hence 

sup \C^{{T,iO,v) \ < sup {J\v\^~^ + J>|"-i) < Ji/^^i-^/^ 

implying that ()2.2()|1 is fulfilled with {Xa, Ha) = («a, 1 — «a)- The corresponding Sobolev expo- 
nent, s < PaC(a/{'ic(a + 2), is then given by, 

1 

S < Si'.-- 



2i + 3 

Finally, for intermediate n's, i < n < 2i, we interpolate the previous two co'£-bounds (which 
are valid for all n's), 

for some ( G [0, 1], which we choose as ( := (n/i) — 1, so that ()2.19|) holds with a = (1 — + 
C/n,pa = (1 - 0/i + 2C/n and ^TMi holds with (A,/i) = (a, 1 - a). This then yields the 
Sobolev-regularity exponent, s < Pa /{3a + 2), 

n+{2e-n)C _ -I, i<n<2i. 

3n + 3{i - n)C + 2nt ^ ^ ' ' 

An additional bootstrap argument improves this Sobolev exponent, s < Pa /{2a + 1), and we 
summarize the three different cases in 

Corollary 4.3. The convection- diffusion equation 

(4-15) |p(^'^) + ^{^P''''(^'-)}-^{;^IA^'^)|p(^'^)}=0' Poe[-M,M], 
has a regularizing effect, po ^ L'^i^x) ^ p{t > e,-) G Wi^^{'M.rc), of order s < se^n given by 



n + {2i-n)Q^n , , If,-. ^ (, ^ ^ 

2n + 2{t — n)Q^n + 




n < 

i < 5 

n>2i 



We note that when n < i, then Qb C Qa and ()4.15|1 is dominated by degenerate diffusion with a 
regularizing effect of order sg^n = Sn = 2/{n + 2). Thus, we recover the same order of regularity 
we met with the 'purely diffusive' porous medium equation (j4.1(Jj) . If n > 2i, however, then 
fla C Qb and it is the hyperbolic part which dominates diffusion, driving the overall regularizing 
effect of (j4.15p with order si^n = = l/(£ + 2); we recover regularity with the same order which 
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we met with the 'purely convective' hyper-Burgers' equation ()3.7|1 . Finally, in the intermediate 
'mixed cases', £ < n < 2£, we find a regularity of a (non-optimal) order 

- ni+{2i-n)in-i) 



We turn to the mult i- dimensional case (j4.1|) . The regularizing effect is determined by the 
size of the set 



1, J > 1, 6<1. 



Assume that the degenerate parabolic part of the equation has a full-rank, so that the smallest 
eigenvalue of b(t>), A(f ) = X^ih^v)) satisfies 

(4.16) I l^b (5) U and sup sup | (b'(i;)^, ^ U ^^,{6) := \v < Xd{h{v)) < . 

In this case, uJc{J; S) ^ {6/ .P)"^^ which yields a gain of VF^^'^^-regularity of order s < 2ab/(2ab + 
1). If in addition, the hyperbolic part of the equation has a non- degeneracy of order a^i namely 

sup |f]a((r,0,'^)U '5°'^ and sup sup |a'(t;) • ^| < 5^-"% 

r2 + |C|2 = l 151 = 1 vm^{5) 

where ila((r, ^); 5) := |i; : |r + a(t>) ■ < then we can argue along the lines of corollary 

14. HI to conclude that there is an overall PF;^'j!-regularity of order dictated by the relative size of 
2ab/(2ab + 1) and aa/(2a;a + 1). As an example, we have the following. 

Corollary 4.4. Consider the two-dimensional convection- diffusion equation 
(4.17) 



s < < 



j,k=l ■' 

with non- degenerate diffusion, B'{v) > |f I". Then, its renormalized kinetic solution admits a 
W[^^ -regularizing effect, po ^ p{t, ■) G W[^^{M.D, of order 

/I 1 \ 

s^m:=niin , , if n > 2max(£,m) and i ^ m, 

^ V£ + 2' m + 2/' - ^ ^ ^ ' 

Sn '■= , ifn<min(i,m) or i = m, 

n + 2 

S£,m,n e [se^rn,Sn], if min(£, TTi) < u < 2max(£,m). 

Finally, we close this section with a third example of a fully-degenerate equation 

In this case, there is a stronger, rank-one parabolic degeneracy with no regularizing effect from 
the purely diffusion part, since (b(f)^,^) = 0, V^i — ^2 = 0, and no regularizing effect from 
the purely convection part where a(f ) • ^ = 0, V^i +^2 = 0. Nevertheless, the combined 
convection- diffusion does have a regularizing effect as demonstrated in 
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Corollary 4.5. Consider the two-dimensional convection- diffusion equation 

For n > 2i it admits a regularizing effect, po ^ L°° (1 L-^ ^ pif^ ) ^ ^zo'c(^^); of order 
s < 6/(2 + 2n-£). 

As before, the Sobolev exponent computed here is not necessarily sharp, and as a result no 
gain of regularity is stated for n < 2i. 

For proof, we use the averaging lemma with the usual {p,q) = (2,1), which yields a 
Sobolev-regularity exponent of order s = [3a{2 — /i — A)/(a + 4 — 2yu), where A and /i 
characterize the degeneracy of the symbol associated with ()4.19|) . 

/:((r, zO, v) = JiT + Jv'iiir + 6) + JM^\ii - 61', + \i? = 1- 

We consider first those ^s, |^| = 1 such that 1^21 < 1/5. Here we have (say) + > 1/10, 
so that VLc{J] 5) G Vta, = {v : \v 1^ < 5/ J}, and (jTTTHl holds with {aa,(3a) = {l/i, 1) : moreover, 
the growth of sought in (j2.2U|) is bounded by, 

sup \C,,{(r,tO,v)\<j(-) +JH-?) ^^V, 

where 

A := l/^ and /i := 1 - l/£, if n > 2£, 

A := 2 - (n - l)/£ and p := (n - l)/£, if n < 2£. 

In particular, if n < 2£ then the Sobolev exponent vanishes since 2 — p — A = 0, and we 
cannot deduce any regularizing effect in this case. If n > 2£, however, we compute the Sobolev 
exponent as before, s = = l/(2£ + 3). Next, we consider the case when — ^2! > 1/5, 
so that nc{J]6) C ^]b = {^^ : l^^l" ^ and (ITTT^ holds with {ab,(3b) = {l/n,2). As for 

^TM . we have 

/ ^ \ (^-1)/" / ^ \ {n-l)/n 

sup \C^{(r,^0,v)\<J(-) + JH-^) ^ J''^', 

where 

A = 1/ra and p = 1 — 1/n, if n < £, 

A = 1/2- (£- l)/n and p = (£ - l)/n, iin>i. 

In particular, if n > 2£ we compute in this case a smaller Sobolev exponent s = 3/(3 — 2£+4n) < 
se. The regularity result follows from the bootstrap argument we discussed earlier which yields 
the final Sobolev exponent s = 2(3a{2 — p — A)/(a + 2 — p). □ 

5. Nonlinear degenerate elliptic equations 
We consider the nonlinear, possibly degenerate elliptic equation 

(5-1) -Y.d^BM^)) = ^^P^'^))^ i^^W, b(-):=B'(-)>0, 

augmented with proper boundary conditions along the C^'^-boundary dV. We assume that the 
nonlinear source term, S'(p), is further restricted so that blow-up is avoided. 
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We begin with formal manipulations, multiplying (j^.lj) against ri'{p), and 'differentiating by 
parts' to find 

=1 i,k=l ■' 



j,k=l 



k 



-trace(v. ® V. B''(p)) + V'(p)(b(p)V.p, V.p) - V(p)5(p). 



We arrive at the entropy inequalities associated with (j^.lj) . stating that sufficiently smooth 
solutions of satisfy, for all convex r^'s, 

(5.2) -trace(v, ® V, B''(p)) < v'ip)S{p), B\p) = [b],{p) := j\'{s)h,u{s)ds]''^ 

So far, we have not specified the notion of solutions for ()5.1|) since it seems that relatively little 
is known about a general stability theory for degenerate equations such as ()5.1|) . The difficulty 
lies with the type of degeneracy which does not lend itself to standard elliptic regularity theory, 
because the BjkS degenerate dependence on p, nor does it admit the regularity theory for 
viscosity solutions, e.g., |ClL92j . |CC95j . because of their degenerate dependence on p rather 
than V xP- We refer to the works of Guan |Gu97j . |Gu02j who shows that in certain cases, one 



is able to "lift" a C^'^-regularity of p into a statement of C°°-regularity. Using the existence 
of smooth viscosity solutions in the uniformly elliptic case where b(f ) > A > 0, ()5.2j) could be 
then justified by the "vanishing viscosity limit" , forming a family of regularized solutions, p^ 
associated with b'^(s) := b(s) + Xldxd and letting A [ 0+. Next comes the kinetic formulation 
of ()5.H) which takes the form 

(5.3) -Vj ■ h{v)VxX,{v) + S{v)-^x,{v) = ^m{x,v) in V'{T x R,), 

for some nonnegative m & Ai^ which measures "entropy production". Indeed, for an arbitrary 
convex "entropy", rj, the moments of ()5.3|) yield 

> — /" ri"(v)m(x,v)dv = [ ■r]'(v)-^m(x,v)dv = 
. . ov 



/d f d 

V(^)&ifc(^) g^,g^^ Xp(^)^^ + J i{v)S{v)—Xp{v)dv = 

J,K=1 ^ 

= - E ^^^Up^ - I {v'{v)S{v)yxpiv)dv = -trace(v. ® V. B''(p)) - V(p)S(p). 

j,k=l 

Thus, the kinetic formulation ()5.3|1 is the dual statement for the entropy inequalities ()5.2|1 . We 
postulate that p is a kinetic solution of ()5.H) if the corresponding distribution function xp{x){v) 
satisfies ()5.3j) . and we address the regularizing effect of such kinetic solutions. 

To use the averaging lemma, we first extend ()5.3|1 over the full x Rt,-space. Let ip be 
C^(-R+)-cut-off function, ipi^s) = 1 for s > e, and let ({x) denote the smoothed distance 
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function to the boundary, ({x) = ip{dist{x,dV)), then f{x,v) := Xpix){v)({x) satisfies, in 



j,k=l -I 



d d d 

—C{x)m{x,v) + ^ -Q-{bjk{v)CxAx)xp{v)) + -Q—{hjk{v)Cx,{x)Xp{v)) + 
d 

S{v)C{x)—Xp{y) - ^ hjk{v)Cx,x^{.x)Xp{,v) 
j,k=i 

d 

(5.4) =: g:^9i{x, v) + A1g2{x, v) + g^{x, v) + gi{x, v). 

Assume that ()5.1|) is non-degenerate in the sense that there exists an a G (0, 1) such that 

(5.5) VL^{5) < (5° and sup sup \{h\v)i,i)\ < d^-"" , VL^{5) := {v e I : {h{v)i,i) < d] . 

We examine the contribution of each of the four terms on the right of ()5.4j) to the overall 
Vl^*' ^-regularity of /, appealing to the different averaging lemmata term-by-term. The first 
term on the right involves the bounded measure gi = (m; averaging lemma ESI with the usual 
(P) q) = (2, 1) then yields that the corresponding average has a VF^^^J ^-regularity of order 
si < 29i, 9i = a/ {3a + 2). The second term on the RHS of ()5.4|) involves the gradient 
of the uniformly bounded term g2 = '^bjk{v)(xk{^)Xp{'v)i here we can use averaging lemma 
12.11 with ?7 = g = 1 to conclude that the corresponding average has VF^*^' -regularity of 
order S2 < O2 = a/ {a + 2) (in fact, with q = 2 one concludes a better 1^*2,2 regularity of order 

52 < a/2). The remaining two terms on the RHS of ()5.4j) yield smoother averages and therefore 
they do not affect the overall regularity dictated by the first two. Indeed, d^Xpi'^) and hence 
g3{x,v) = S{v)({x)dvXp{'v) is a bounded measure and averaging lemma ITT] with 77 = = 
implies that the corresponding average /g belongs to the smaller Sobolev space, Wi^^ of order 

53 < 26^3, 6^3 = a/{a + 2). Finally, the last term on the right of ()5.4|) consists of the bounded 
sum, g4 = — J2bjk{ij)CxjXk{x)xp] with r] = N = and q = 2, the corresponding average /4 has 
a W[^f^ regularity of order S4 < 26*4, 64 = a/2. 

Next, we iterate the bootstrap argument we mentioned earlier in the context of hyperbolic 
conservation laws. The first W^^^'^-bound together with the L°°-bound of / imply a PF^'^^'^-bound 
with (Ti = si/2, which in turn yields the improved regularity of fi G PF^q^ , s < (1 — 9i)ai + 29i. 
Thus, for the first term we can iterate the improved regularity, si (1 — 9i)si/2 + 29i, 
converging to the same fixed point we had in the parabolic case before, si < 2Q;/(2a + 1). 
The second term requires a more careful treatment: as we iterate the improved regularity of 
Xp ^ ^loc^ '^^^ express the term on the right of ()5.4|) as A'^''g2 with rjs := I — s and with g2 
standing for the sum of L^-bounded terms, g2 = A'^ '^bjk{v)(xkXp{'>^)- Consequently, averaging 
lemma ini vields the fixed point iterations S2 ^ (l — 6'2)s2/2 + (2— 77^2)^2 with limiting regularity 
of order S2 < a. The remaining two terms are smoother and do not affect the overall regularity: 
a similar argument for the third term yields the fixed point iterations, S3 (1 — ^3)53/2 + 2^3 
with a fixed point S3 < 2a/{a + l), while the fourth term remains in the smaller Sobolev space 
VF"'^. We summarize with the following statement. 
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Corollary 5.1. Let p G L°° be a kinetic solution of the nonlinear elliptic equation 45. ij) and 
assume the non- degeneracy condition / I5.5)) holds. Then we have the interior regularity estimate 
for all D CT, 
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